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Examples of Systems Demanding Predictive Models of Effects of Gravity on 

Two-Phase Flow and Heat Transfer 
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NASA-Supported Facilities at Boiling & Two-Phase Flow Laboratory 
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PURDUE Overriding Objectives of FBCE 

university o 


The proposed research aims to develop an integrated two-phase flow 
boiling/condensation facility for the International Space Station (ISS) to serve as 
primary platform for obtaining two-phase flow and heat transfer data in 
microgravity. 


Overriding objectives are to: 

1. Obtain flow boiling database in long-duration microgravity environment 

2. Obtain flow condensation database in long-duration microgravity 
environment 

3. Develop experimentally validated, mechanistic model for microgravity flow 
boiling critical heat flux (CHF) and dimensionless criteria to predict minimum 
flow velocity required to ensure gravity-independent CHF 

4. Develop experimentally validated, mechanistic model for microgravity 
annular condensation and dimensionless criteria to predict minimum flow 
velocity required to ensure gravity-independent annular condensation; also 
develop correlations for other condensation regimes in microgravity 
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Purdue 

UNIVERSITY 


Layout of FBCE 



Consists of: 

■ nPFH sub-loop 

■ Water sub-loop 

Contains three test modules: 

■ Flow Boiling Module (FBM) 

■ Condensation Module CM-HT for heat 
transfer measurements 

■ Condensation Module CM-FV for flow 
visualization 
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PURDUE I Flow Boiling Facility 

UNIVERSITY 
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Flow Boiling Module 
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Heated Wall Design 
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Horizontal Flow Boiling - Heated Wall Con figurations 

UNIVERSITY 
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Horizontal Flow Boiling - Flow Visualization 


Slightly subcooled flow at low velocity 


Top Wall 
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Horizontal Flow Boiling - Experimental CHF Results 
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Separated Flow Model for Double-Sided Heating 


Mass Balance: 
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UNIVERSITY 


Interfacial Lift-off CHF Model 


Use separated flow model to determine axial 
variations of: 
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PURDUE I CHF Model Predictions 

UNIVERSITY 


One g e Horizontal Flow 


Microgravity & 
One g e Vertical Upflow 




0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 1.6 1.8 2.0 


U [m/s] 


U [m/s] 



ASGSR 2015 


Boiling and Two-Phase Flow Laboratory (BTPFL) November 2015 15 





Purdue 

UNIVERSITY 


Single-sided versus Double-sided Heating in Earth Gravity 
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CHF Predictions for Single-sided versus Double-sided Heating in Earth Gravity 

UNIVERSITY 
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Data Sharing Plans 
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Organizational Structure of FBM Database 


Folder Name: 
Summer 2015 
FBM Testing 



ASGSR 2015 


Boiling and Two-Phase Flow Laboratory (BTPFL) November 2015 


20 





Purdue 

UNIVERSITY 


Data Folders 


Contain four filetypes: 

■ Text Files containing raw sensor data output by data acquisition system 

■ Matlab Scripts for processing raw data 

■ Excel Spreadsheets containing all relevant parameters (e.g., pressure drop, heat 
transfer coefficient, CHF) output by processing script 

■ Image Files for flow visualization 

With subfolders used to group data by operating conditions 



T&B, 2 Phase, 
3psig Accum 
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• • • • 
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found in 
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Documents” 
folder 
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Computational 

Modeling 
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PURDUE ■ ANSYS Fluent Modeling 

UNIVERSITY 
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Numerical Approaches to Modeling Two-Phase Systems 

UNIVERSITY * * 


• Lagrangian 

- Smoothed-Particle Hydrodynamics (SPH) Method: Gingold & Monaghan 

(1977), Lucy (1977) 

- Multiphase Particle-in-Cell (MP-PIC) Method: Harlow (1955) 

• Eulerian 

- Level-Set Method (LSM): Osher & Sethian (1955) 

- Volume-Of-Fluid (VOF) Method: Hirt & Nichols (1981) 

- Coupled Level-Set/Volume-Of-Fluid (CLSVOF) Method: Sussman & Puckett 
(2000), Enright et al. (2002), Tomar et al. (2005) 

• Eulerian-Lagrangian 

- Front Tracking Method: Unverdi & Tryggvason(1992), Tryggvason etal. (2001) 
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Phase Change Models 


Rankine-Hugoniot jump condition 

Mao (2009). Michita & Thome( 2010), Sun etal. (2012) 
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Phase Change Models 


Rankine-Hugoniot jump condition 

Mao (2009). Michita & Thome (2010), Sun etal. (2012) ... 

Pros: Cons: 

1. Ease of implementation 1 


2 . 


Allows for phase change only along 
interface 

Cannot maintain saturation 
temperature 


Schrage Model (1953) 

Kartuzova & Kassemi (2011), Magnini etal. (2013)... 

Pros: Cons: 

1. Successfully used for 1. Requires use of empirical coefficient y 

evaporating & condensing films 2 A||ows for phase change on|y a|ong 

interface 


Lee Model (1980) 



Wu etal. (2007). Yang etal.{ 2008), Fang etal. (2010) ... 

Pros: Cons: 

1. Ease of implementation 1. Not applicable for subcooled boiling 

2. Successfully used for 2. Requires use of empirical coefficient r ; 

condensation processes 


ASGSR2015 Boiling and Two-Phase Flow Laboratory (BTPFL) November 2015 26 





Purdue 

UNIVERSITY 


Computational Domain , Governing Equations and Boundary Conditions 


Computational Domain 



Governing Equations 

■ Continuity Equations 

- Liquid Phase: 

^(«/ p /) + Vx (“/ p / s /) = 0 

- Vapor Phase: *- ee m °d e l (2008) 

^(“,p ! ) +Vx (“,p«" s ) = 0 s r- s f= r 


* g g j 1 
sat 


Momentum Equation 

— (ptt) + Vx(p««) = -V.P + Vx p^V« + V« r j 


+ pg + F 


Energy Equation 

^(p£)+vx(s(p£ + i>))=vx(* !(r vr)- l g<^a=/, /i; s / 


Boundary Condition 

■ Axisymmetric centerline 

■ k-u SST turbulence model 

■ Inlet uniform velocity from experimental data 

■ Wall heat flux from experimental data 
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Void Fraction Results: Climbing Film Regime 
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Predictions of Average Heat Transfer Coefficient 
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Future Computational 

Modeling 
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Purdue 

UNIVERSITY 


Current Work 


Current Work 

■ Working with Prof. Carlo Scalo’s 
group at Purdue University to 
develop a 2-D code using best 
available computational techniques 

- Performing comparisons with Fluent 
to quantify in-house solver 
performance 
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Future Work 
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Use proposed 2-D code to run select cases which can be represented 
reasonably well by 2-D domains (e.g., axi-symmetric flow condensation, 
flow) 

Scale 2-D code up to 3-D, highly parallelized solver, which will include 
turbulence effects, to be run on Purdue supercomputing cluster 

Begin comparing data for transient cases with 3-D geometry to prior 
experimental studies 
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